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Применение метода сечений для контроля 
формы поверхности пятна излучения 
в реальном времени 


В статье рассматривается проблема контроля формы поверхности пятна излучения в реальном времени. 
Предлагается использование метода сечений для формирования выборки и последующей обработки 
дальнейших изображений. 


Введение 


В настоящее время в полиграфии, лазерной обработке материалов, локации, опти- 
ческой связи и других областях техники ощущается необходимость более широкого 
внедрения оптико-электронных систем с автоматической корректировкой искажений 
формируемого светового излучения. Причинами этих искажений могут быть дестабили- 
зирующее воздействие механических или климатических факторов, нестабильность 
характеристик источника излучения, возмущения в оптическом тракте, разъюстировка 
оптических элементов и т.п. Обеспечение приемлемого качества коррекции требует 
непрерывного динамического контроля характеристик светового излучения, например, 
пространственного распределения его интенсивности, в том числе оценки отклонения 
указанного распределения от исходного или эталонного распределения. 

Целью данной работы является разработка метода выборки и последующей 
обработки формы поверхности пятна излучения в реальном времени. 


Приближенная классифицирующая функция 


Традиционный путь решения данной задачи включает формирование изображения 
пятна излучения В\(х, у) на светочувствительной поверхности фотоприемника с после- 
дующим преобразованием в сигнал Ох, у), амплитуда которого в каждой точке разло- 
жения с координатами (х, у) соответствует интенсивности в В (х, у), т.е. Ибх у) <> В(х, у). 
Далее сигнал С(х, у) сравнивают с эталонным сигналом И/(Х, у) для всех точек разложения 
сигналов. Сигналы ((х, у) и Их, у) представляют собой некоторые поверхности, которые 
могуг отличаться типом, коэффициентом относительного масштаба, вектором относи- 
тельной координатной привязки и углами относительных поворотов в трехмерном 
пространстве. Поэтому сравнение этих поверхностей должно производиться с уче- 
том всех возможных ситуаций, что требует огромного объема вычислений и затруд- 
нительно для осуществления в реальном времени. 

В практически важных случаях требуемое сравнение поверхностей И(х, у) и 
И/Хх, у) может быть реализовано с помощью излагаемого ниже нового метода сечений. 
Этот метод включает следующие операции: 
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1) нахождение максимальных амплитуд сигналов Ц(х, у) и Их, у) (рис. 1): 
О „ = шах О(х,у)=Й,, 


тах 


Й к = ах Й(х,у)=Й,; 


тах 
2) нахождение значений площадей 5,.;, и 5,., сечений на уровнях 0,75 Й, и 


0,5 И», для сигнала ((х, у) соответственно, бу-;„ И 5%°5„ на уровнях 0,75 й» и 0,5 Й„, 


для сигнала И/х, у) соответственно; 
3) расчет приближенных значений коэффициентов формы 7,, и 7, „ для сигналов 


ОС, у) и И’, у) соответственно; 


пи, ЕК = Зови / Зобь (1) 
пи бы ЕК» = бон (Зоб › (2) 
где - — знак приближенного значения, 
Пе (3) 
И И Е (4) 


И и и И, - суммарное значение амплитуд сигналов И(х, у) и Их, у), меньших 
уровней 0,5 Й, и 0,5 й,„, соответственно. 

Покажем справедливость 7, = Г, на примере текущего сигнала ((х, у) следующим 
образом. 

Объем ра фигуры выше сечения поверхности ((х, у) по уровню 0,5 й, может 


быть найден с помощью одной из известных формул приближенного вычисления 
интегралов для равноотстоящих узловых точек, например, формулы Симпсона [1, с. 507]: 


Е 45 5 ес 5 0.5 ]/12 ы Е (4. + 1) /12 . (5) 


Г и & й, [5(( 


тах 


Так как 5(И 


амплитуды, при этом чаще всего 5, #0; 


)=5„- площадь сечения поверхности на уровне максимальной 


тах 


А у/й, у\ 
и, 


Рисунок 1 — Сигналы эталонного и текущего изображения пятна излучения 
4) сравнение коэффициентов формы 7. ии. „ (вместо поэлементного сравнения 


поверхностей). 
Используя (5), легко получить 


Г, и = Г, и нЕ ЙО би /2 ы Я оби (47. + 7)/12 2 (6) 
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откуда следует 
и = (47,,+7)/12. (7) 
Линейная зависимость коэффициентов 7, и г, позволяет использовать коэффици- 


ент г, как характеристику формы соответствующей поверхности. 
Возможности классификации типов поверхностей посредством коэффициента 
7, отражены в таблице. Разные значения коэффициента 7, соответствуют разным 


типам поверхностей приведенных в ней примеров фигур. 


Таблица 1 — Профиль продольного сечения 


Значения коэффициента Значение 
Рисунок продольного погрешности 
сечения и. 7, и, бр 
1. Призма, цилиндр . 
А и 
т Е, 1 1/2 = 0.92 1 1/12 = 0.083 
2. Пирамида, конус 
и 
1,0 
_ 14 =0.25 | 2/3=0.67  2/3=0.67 0 
0,5 "В 
Хх 
7/12 =0.58 7/9 =0.78 7/9 = 0.78 0 
1 О: 
4. Полушар 
с основанием 
2 Тв 
ть 34=0.75  5/6=0.83 | 5650.83 0 
1,0 
> 
0 1 Х/К 
5. Гауссоида 
(//2л)ехрЕ(х? + у2)/2] 
А д 
и | 1/3) 12 | 66/12 = 0.72212 ш=0.721 =0 
ТЕХ = 0.415 
0,5 
0 ” 
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Важным достоинством коэффициента г, для некоторых типов поверхностей явля- 
ется его независимость от коэффициента масштаба, сдвига и ориентации соответству- 
ющей поверхности. Кроме того, вычисление коэффициента г, ‚ так как число разрядов 


кода значения 5, всегда существенно меньше числа разрядов кода значения Г, .. 


Коэффициент 47, характеризует обобщенную выпуклость поверхности: при 
4х. >1 — поверхность выпукла, при 4и, <1 — поверхность плоская в обобщенном 
смысле, при 4и. =1 — поверхность линейна в обобщенном смысле. 


Важно отметить, что поверхность можно охарактеризовать площадью эффектив- 
ного сечения 5.: 


й, = 2(7% 5 — 50-й -0,5) /й = 2(5,57й-5,.-й.0,5)/й =25,. (7, -0,5), (8) 
откуда 
(45 5 +5%5)/6 при 5, =0, 


Я =$ = 
у ы (5, +450 5 +505)/6 при 5, #0, 


(9) 


где 5, — площадь сечения сигнала на его максимальном уровне. 


Метод сечений [2] можно использовать для выборки сигналов лазерных пучков, а 
в качестве параметра выборки использовать коэффициент формы, при этом данная 
методика будет использоваться для нахождения центров лазерных пучков. 

Алгоритм заключается в следующем: 

1. На изображении находится точка с максимальной яркостью. 

2. От найденного значения яркости на изображении выбираются 5 значений 
яркости, каждое на единицу меньшее от предыдущего. Крайние точки каждого из 6 
диапазонов образуют краевые линии. В результате получается 6 краевых линий. 


Г 
«ЕЕ 
ИЯ 


1% 
> 
я 


а) 6) 
Рисунок 2 — Форма пятна лазерного пучка 
(а — фото с камеры, 6 — схематический рисунок) 


3. Для каждого краевого контура определяется предварительный центр, после 
чего вычисляется усредненный центр для координат Х и У. 

Координаты (х,у) энергетического центра двумерного сигнала с удельной нели- 
нейной плотностью и(Кх,у)), соответствующей данному значению /(х,у), выражаются 
в следующем виде: 


МЬ-—1М- М-—1 №М- 


1 №-—1 1 №М-—1 №М-—1 №М-—1 
ли Уж еу))-х, у= =. Уж бу), М= >>» 06). 
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4. Полученные краевые линии разбиваются на 4 сегмента. В качестве центральной 
точки используются координаты предварительного центра (усредненного центра). 

5. Определяются «хорошие» и «плохие» изображения по факту попадания в 
рассчитанный туннель вычисленного коэффициента формы. Если коэффициент формы 
текущего изображения попадает в туннель, то последнее считается «хорошим», в 
противном случае — «плохим». Границы туннеля вычисляются в обучаемой выборке 
изображений путем последовательного выбора изображений и соответствующих 
коэффициентов формы с минимальным разбросом координат предварительных центров. 

6. Вычисляются координаты энергетических центров только для «хороших» 
изображений. 


Результаты экспериментальных исследований 


Исследовано 15 трасс лазерных пучков [3] и получены следующие результаты 
(показано для 2 трасс): 

1. Границы туннеля а := 0.777 Ъ := 0.78. 

Значение разницы между максимальным и минимальным значением координат 
1.147. 


0 200 400 600 800 


Рисунок 3 — График распределения координат центров лазерных пучков 
(а — значение Х, б - значение У), где Х, У - координаты всей трассы лазерных 
пучков соответственно, ХМ, УМ -— значение координат центров пятен после выборки 
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2. Границы туннеля а := 0.751 Ъ := 0.757 
Значение разницы между максимальным и минимальным значением координат 2.18. 


0 200 400 600 800 


0 200 400 600 800 


Рисунок 4 — График распределения координат центров лазерных пучков 
(а — значение Х, б — значение У), где Х, У - координаты всей трассы лазерных 
пучков соответственно, ХМ, УМ -— значение координат центров пятен после выборки 


Выводы 


Эксперименты показывают, что для изображений пятен лазерных пучков в резуль- 
тате действия разнообразных дестабилизирующих факторов невозможно точно изме- 
рять координаты их энергетических центров, однако точность можно существенно 
повысить, используя вычисление коэффициента формы изображений с последующей 
их классификацией на «хорошие» и «плохие». Сравнительный анализ показывает, что 
предложенные методы позволяют измерять координаты центра лазерных изображений 
с точностью определения точки привязки не менее чем 1 — 2 элемента разложения, что 
по точности превышают известные, например, на основе определения центра масс с 
помощью моментных признаков, в среднем в 1,5 раза. 

В соответствии с изложенным методом сечений сравнение поверхностей сводится 
к сравнению их коэффициентов формы и не требует поэлементного сравнения с 
учетом всех случаев отличия их типов, коэффициента масштаба, относительного 
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сдвига и поворота в пространстве. Важно отметить, что равенство коэффициентов 
формы поверхностей, в общем случае, позволяет отнести их к одному и тому же 
обобщенному типу поверхности или аппроксимировать этим типом поверхности. 

Очевидно, что метод сечений может быть легко распространен на случай увели- 
чения числа равноотстоящих сечений поверхности или изменения значений уровней 
сечений, изменения весовых коэффициентов площадей сечений. 

При этом в качестве порождающей основы также целесообразно использовать 
подходящую формулу приближенного вычисления интегралов [1, с. 507]. 

Рассмотренный метод сечений перспективен также для применения в задачах 
классификации и архивации изображений в реальном времени. Важным достоинст- 
вом данного метода является простота осуществления как программными, так и 
аппаратными средствами. 
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Застосування методу перетинв для контролю форми поверхн! плями випромнювання 

в реальному час! 

У статт! розглядаеться проблема контролю форми поверхн! плями випром1нювання у реальному час1. 
Пропонуеться використання методу перетивв для формування вибрки 1 наступно! обробки подальших 
зображень. 
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